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o 

Cn ' We study the role of electron-electron interactions for spin polarized transport using the Boltz- 

mann equation and derive a set of coupled transport equations. For spin polarized transport the 

^ ' electron-electron interactions are important, because they tend to equilibrate the momentum of the 

two spin species. This "spin drag" effect enhances the resistivity of the system. The enhancement 

>0 ■ is stronger the lower the dimension and should be measurable in for example a two dimensional 

electron gas with ferromagnetic contacts. We also include spin flip scattering which has two effects: 

it equilibrates the spin density imbalance and, provided it has a non s-wave component, also the 

jrt ' current imbalance. 
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I. INTRODUCTION 



Recent advances in fabrication of ferromagnetic-semiconductor hctcrostructures |l| and the observation of spin 
C^ injection into semiconductors [^ have lead to new interest in transport properties of spin polarized systems. There 
has been considerable work done in the field of metallic magnetic multilayers, which has been analyzed in terms of 
I ' transport equations with spin dependent distribution functions |||,Q|. These works based their analysis on diffusion 
' O ' transport equations. The justification of using these equations was given by Valet and Fert Iq], who derived a spin 
^ diffusion equation from the Boltzmann equation in the limit where the spin scattering length is much longer than 

the momentum relaxation length. Recently, the transport equation were utilized to analyze the feasibility of spin 
injection into semiconductors with the result that the crucial parameter is the conductivity mismatch between the 
semiconductor and the ferromagnet B and to study theoretically spin-polarized transport in in inhomogeneous doped 
semiconductors 0. 

None of these approaches take the electron-electron (e-e) scattering into account. Clearly the e-e interactions play 

^|J ■ a different role than in usual spin degenerate transport, where the e-e interaction does not provide a mechanism for 

,_^ momentum relaxation and hence has only indirect consequences for transport coefficients. In spin polarized transport 

f — I the two spin species have different drift velocities and e-e interactions are instrumental in equilibrating this difference. 

^^ ' This leads to a spin drag effect where the spins carrying the larger current will drag along the spins carrying the 

smaller current. This drag effect has been considered recently by D'Amico and Vignale |§,H| in three dimensions using 

linear response theory. They found that the the spin drag resistivity was appreciable and can at elevated temperatures 

become fractions of the usual resistivity of the metal. In two dimensions the effect of e-e interactions on spin diffusion 

has been considered theoretically by Takahashi et al. |10| , |ll| . Using a quantum kinetic equation approach, previously 

utilized in ■^He-'*He solutions |l3], they studied the spin diffusion coefffcients in two dimensional electron gases. In 

i-rt ■ order to study this spin diffusion they used variational functions but did not include spin relaxation scattering. 

^ ' In this paper, we use the Boltzmann equation to study spin dependent transport and spin diffusion. We restrict our 

O ■ selves to the study of colinear magnetization and our goal is to derive a set of transport equation in the semiclassical 

O limit. For this purpose the Boltzmann equation is adequate. For the non-colinear case, where phenomena like damped 

J> , transverse spin modes cannot occur, one must go beyond the present approach, see e.g. [|l^,0 and references therein. 

We include impurity scattering, both spin independent and spin flip scattering, as well as e-e scattering. We show 

that a shifted Fermi-Dirac distribution (SFD), is a valid solution at low temperatures, T <^ Tp and without spin 

jrt ■ flip scattering. This is also the case for weak e-e scattering, where the problem in absence of spin flip reduces to the 

" " ' ordinary Coulomb drag situation [f3UlJ] . 

We then go on to discuss the general case at higher temperatures, general interaction strength and flnite intrinsic 
spin flip scattering. Using a SFD ansatz, we find for an isotropic system the following macroscopic transport equations 

/ e dn^\ 



V-J.= (^-^^j(M.-M-.), (la) 

Va2s = — J. + (— ^ + -^ ) (Js-asJ-s)- (lb) 

1 



Here Js is the current carried by electrons with spin s, jig is the local spin dependent electro chemical potential, as 
is the conductivity of the spin s electron gas, t°£ is a spin lifetime due to intrinsic spin flip processes, asi.s is a spin 
current conversion conductivit y a rising from the angle dependence of spin flip scattering, and as = n^g/rfi_g is the 
relative spin density, see Eqs. (pq) , (|3S|) for definitions. Finally ajj is the spin drag conductivity given by 



{^d)-' = 



de2 



HsU^s 



(27r)'' 



dujq'^\e(j){q)\ 



2 lmxs(g,tx;)lmx-s(g,^) 
kBTsmh'^(nuj/2kBT) 



where Im Xs is the polarization function 



/dk 
j^ [/°(|k + q|) - fs{k)] S (£|k+q|,. - Sks - tiw) 



(2) 



(3) 



The formula (0) is well-known from Coulomb drag |jT3]. At low temperatures it goes (an)^^ oc T^ in two and three 
dimensions, while in one dimension it is proportional to T, see e.g. |l5[ . 

The first transport equation (y_aj) is the continuity equation, which expresses the conservation of current in the 
presence of spin flip processes. The second equation (lb) is a generalized Ohm's law. The first term on the right 
hand side is the usual Ohm's low, while the second term describes that a momentum imbalance between the two 
spin directions gives rise to an additional resistance if there is a mechanism for conversion of spin current. There 
are two such processes possible. This first one is the spin drag effect mentioned above, where e-e scattering makes 
transfer of momentum possible. The second one is due elastic spin flip scattering on for example magnetic impurities, 
which can convert a current with one spin polarization to a current of the opposite polarization, if the spin flip matrix 
element has an angular dependence. For example if the spin flip predominately scatters forward it means that a 
spin flip scattering is accompanied with a transfer of momentum. In contrast if the spin flip scattering is purely 
s-wave scattering the momentum transfer between the spin channels is on average equal to zero. This can be seen 
mathematically from the expression for Tgf in Eq. (35c). The derivation of the these two terms is the main result of 
the present paper. 

Two consequences of the spin current relaxation terms can immediately be read off. Firstly, they give rise to an 
increased resistivity in the case where the current is spin polarized. For example, taking J^ = the effective resistivity 



for electrons with spin f becomes (l/cr| + l/aoag + 1/cr, 



Eqs. (|l|) we get a diffusion equation for the electro chemical potential difference 



and hence an enhanced resistivity. Secondly, from 



where 
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(5) 



This shows that the intrinsic spin relaxation length is decreased by the spin drag and angle dependent spin flip effects. 
Similarly, we get that the following weighted sum of electro chemical potentials must vanish. 



V2(c_sAis + Csfl-s) = 0, 



(6) 



where 



9ng 1 
di-J. (Js 



Cst.s 



(yocts 






(7) 



Below we derive Eqs. (|la|) , ([lb|) and estimate the spin drag contributions. For the two-dimensional case, we have 
also performed the integration of Eq. (0) numerically. 



II. BOLTZMANN EQUATION FOR COLINEAR SPIN TRANSPORT 



We base our analysis on the Boltzmann equation for transport through a system with lifted spin degeneracy. We take 
the current to run in the a;-direction and denote the non-equilibrium distribution function /s(k) and the equilibrium 
Fermi-Dirac distribution function /g 



where ^o is the chemical potential and [3 as usual the inverse temperature. The eigenenergies are denoted e^s, where 
s is the spin quantum number and k the quantum number labeling the relevant states crossing the Fermi level. For 
simplicity we assume a parabolic dispersion and we write 

eks^%^+e1, (9) 

where e^ is the band off-set which can be spin dependent if the material is ferromagnetic. 
The linearized Boltzmann equation then reads 

We take the collision integral to include elastic scattering and e-e scattering 



= ifo[/,s](k)+ff.f[/„/-.](k)+i/e-e[/.,/.](k)+He-e[/s,/-s](k), (11) 

where Hq is the scattering from impurities (or quasielastic phonon scattering) giving rise to a momentum relaxation 

r r/k' 
i7o[/.](k) = -j ^^M/0(k,k') [fs{k) - fs{k')] S{e,s - e,,.), (12) 

and where Hgf describes elastic scattering processes that flip the spin 

/r/k' 
^2^Ty,/(k,k') [/,(k) - /_,(k')]<5(£fc, -efc,_,). (13) 

Finally, the e-e scattering is after the linearization given by |16| 



ife-e[/.,/.'](k)--^y^^y^|f/..,(q,£k-ek4 



d\' 



x^(ek + ek' — ek+q — £k'-q) 

x.fo{£ks)fo{ek's') [1 - /o(£|k+q|s))] [1 ~ /o(£|k'-q|s')] 

X [*,(k)-f *,,(k') - *,(k + q) - *,,(k' - q)] , (14) 

where the deviation from equilibrium is expressed in the function '^ through 

/.(r,k) ^/°(fc)+ (-^^) *.(r,k). (15) 

The interaction Uss' is the Coulomb interaction between two electrons with spin s and s' . It can in principle depend 
on the relative direction of the spin if exchange is included. This set of integral equations cannot be solved in general 
and one must either solve them numerically (for example as in Ref. [pT|) or proceed with approximate methods. 

One simplification is however possible from symmetry. Because of the cylindrical symmetry the functions ^s(k) 
only depend on the angle between k and the direction of the current, which we here choose to be in the x-direction. 
Denoting this angle by 9 we have cos 6* = k ■ x/A: and we can write 

^,ir,k)^^,{x,k,9). (16) 

It is convenient to expand the distribution function in harmonics of the angle 9 as 

■^ s{x, k, 9) ^J2 9 s''Hx,k) cos ne, (17) 

n=0 

which we utilize in the next section. 



III. SPIN DRAG WITHOUT SPIN FLIP PROCESSES FOR T «.Tf 

In this section we study the Bohzniann equation in the presence of e-e interaction, but in the absence of spin 
flip processes, i.e. H^f = 0. Furthermore, because a low temperature expansion allows for a solution of the Boltz- 
mann equation, we start by examining this limit and later we discuss the validity of this solution even at elevated 
temperatures. It turns out that the solution in the low temperature regime corresponds to a SFD distribution. 

In the low temperature limit we see from Eq. (|lO|) that the second term on the left hand side (the driving term) 
restricts Sks to lie close to the Fermi level, such that the deviation ^s (k) need only be evaluate at kp ■ This is therefore 
also true for the distribution function in the elastic collision term, Hq. Due to the Pauli principle this will also be the 
case for the ^'s in the e-e collision integrals, which is seen as follows. Using standard tricks (see e.g. [H), we rewrite 
the e-e collision term as 

^e-4/./.](k)^-f£g/(^/(0,|f/..(q,e.-e,,,)P 

^Tsinh^(L/fc.T)'"^^^('''^'")''"^^'('^'-^="^ 

X [vI/,(k)-KvI/,,(k') - vl/,(k + q)-*,,(k' - q)] , (18) 

where 

Imxs(k,q;tj) = tt [/o(e|k+q|J - /o(efcs)] ^ (e|k+q|s - £ks - M ■ (19) 

Now, at low temperatures the factor 1/ sinh restricts the tj-integral to small w of order fcT, and hence e|k+qis in ( |l9| ) 
deviates from ets by an amount of order kT and we expand Im Xs as 

Imxs(k,q;tj) « tt^w f " ^^ j (^ (£|k+q|s - £ks - fiuj) ■ (20) 

From this we conclude that both Sks and et's' (and hence also e|k-i-q|s and £|k'-q|s') are within a shell of order fc^T 
from the Fermi level. To leading order in kT/ep we can therefore neglect the dependence on k and keep only the 
angular dependence of ^g. Therefore in the following we replace 

*,(k)«*,(fcF3,0), (21) 

where kps is the Fermi wave vector for the spin direction s. 

Now we expand the function \1/ in harmonics of the angle 9 as in Eq. (|l^) . Inserting ( p^ ) and (|7|) into the Boltzmann 
equation gives for the left hand side 

Tifcsx / dfQ{eks)\ fd'^s „\ >r-^ Tifc^ cos g / dfo{eks)\ f ^dg^"^^ .^ , 

e£;, =>^ \\cosne—- (22) 



dsks J \ dx J ^-^ m V dsks J \ dx 



and for the right hand side we have two terms. The first one is the spin conserving impurity scattering term which 
becomes 

where we defined the transport times of order n 

1 /■ dk' 

^^= / 7:rT;jl^r(k,k')[l-cos7^0kk']<5(£fe.-£fe-.), (24) 

and where ^kk' is the angle between k and k'. The second term is the one with the e-e scattering. When the expansion 
( p^ ) is inserted into the e-c interaction terms different n do not couple, see for example the derivation in Ref. p7[ . 
The trick is to write for example the angle of k' — q as cos ?T.0k'-q,x = cos n (0k'-q.k + ^k,x) — cos 7i^k'-q,k cos rt^k,x ^ 
sinn0k'-q,k sinn^k.x and note that the sin terms vanish due to symmetry. Therefore we can express the e-e collision 
term as 



H,^,[fs,fs']{k,0) =^cosM) (5i"^j(") +5lrV(")) , (25) 



where J*^"^ corresponds to the first and thnd term in Eq. ( |l^ ) while /f^"^ corresponds to the second and fourth term. 

(n) 

Now a set of equations for the coefficients gs can be extracted by multiplying the Boltzmann equation with 
cosn'9 and integrating over 9, while using that J d9 cos n9 cos n' 6 oc (5„„/. The left hand side of Eq. (E2h is expanded 
in harmonics using that cos 6* cos n6' — ^ [cos6'(n + 1) + cos9{n — 1)]. We find the following set of equations 
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(26a) 



|tiL(,(.), 5p, ,,),.= (,a.|, ,,<.,(,;<«, ,<.,),,.«,«), „6b) 



(9:"-"+rf'«')').= rf"|j+9S"'(j'"'+/"")+9l".'/"" , nSa (26c) 



where 



The solution of these equation is that gs — for n > 2. This is due to the fact that dxgs = which decouples 
the equations for n > 2 from the two first. Eq. (E6a) expresses current conservation within each spin species. If we 
include spin flip scattering in the equation, then the equations couple because dxgs 7^ 0. 

Now we note that setting 17^ = for n > 2 corresponds precisely to a linearized shifted Fermi-Dirac (SFD) 
distribution 

/f™(fc) = /o(e.(k + (5k,)-JM,), (28) 

from which we read off (for (5k in the x-direction) , 

4-2 7 

gf^^5^ls. g^^^cos9^nvJks^g'^P ^ —5ks. (29) 

ra 



From Eqs. (26), one can now determine 5^ and g].. They correspond to the change of the local charge densities and 
to the local currents, respectively. This we will see in the next section, where we use the SFD ansatz to study the 
general case. 

IV. MACROSCOPIC TRANSPORT EQUATIONS 

Above we saw that at low temperature the exact solution of the Boltzmann equation in absence of spin flip 
was a shifted Fermi-Dirac function. The same conclusion applies to the situation of arbitrary temperatures but 
weak e-e scattering, because this limit corresponds to the usual Coulomb drag regime. This is however no longer 
necessarily true at arbitrary e-e scattering, when the temperature increases, or when spin flip processes are included. 
Nevertheless, we shall assume in the following that the SFD distribution is a good approximation for the exact 
distribution function. The argument for doing this is as follows: the e-e interactions will drag the distribution functions 
towards shifted Fermi-Dirac distributions because the inter spin channel e-e collision terms vanish for fs = ff^^, i-e. 
He-e[ff^^ , ff^^] = 0. Since the e-e scattering rate increases as rJL^g « {eF/h){kT/eF)'^, increasing the temperature 
actually helps. Furthermore, since the energy dependence of the elastic scattering is important in determining the 
actual shape of the distribution functions, and because we do not want to into details of this sort, we will use view 
the SFD distribution functions as reasonable parameterizations of the true distribution function. 

Our starting point is thus an ansatz distribution function given by 

fs{k) = Me,s) + [-^^^) Sf^six) ~ (-^1^) nvxSksix). (30) 

Here Sfig corresponds to a change of the local chemical potential and hence also to the local density, while 6ks describes 
a shift of the distribution function in k-space and thus gives a finite drift velocity. Inserting this into the Boltzmann 
equation gives for the left hand side 



LHS = v^ (-^1^BA\ fl. (5^^ _ nvjks) ~ eEj , (31) 

and for the right hand side we have three terms. The spin conserving coUision term becomes 

Holfs] - -hvjk,^ fdMeksy\ (32^ 

Ttr,s V (J^ks J 

where the usual transport time is 

1 f dk' 

= / TT^^siK'i^') [1 - cos^kk'] S{eks - Sk's). (33) 

TtT,s J 1^71'J 

The second term on the right hand side is the spin flip scattering term which becomes 

HsiifsJ^s] = Uvj:Sks {Sfis - Sfi-s) -n- - ^v^ {Sks - Sk^s) — I T^ — ] , (34) 

V Tsi,tv T^f Tsf/ V 0£ks ) 

where the three different spin flip scattering times are given by 

^=y (^M/./(k,k')5(£..-s.,_.), 

"^ / 7:5-Td^-'/(k,k')[l-cos6ikk']<5(efcs ~efe'-s), 

1 f dk' 

— ^ / 777^:7 ^5/ (k,k') cos 6'kk' ^(^fes -£fc'-s)- 
Tsf J (^tt)" 

Finally, the e-e scattering is given by J^s' He-e[fs,fs']- But in accordance with detailed balance the e-e scattering 
between two identical Fermi-Dirac distributions is zero, Hf,-e[fs, Is] ~ 0, and we are left with 



(35a) 
(35b) 
(35c) 



Tj r^ ^ 1 h' f°" dio f dk' f dq . ^ , ,,2 

m J-oo 27r J {2TTf J {2TTf 



X- 



■Imxs(k,q;cj)fmx-s(k',-q;w) 



kT sinh (Huj/kBT) 

xq^[Sks-Sk_s]- (36) 

The final form of the Boltzmann equation is thus 

(31) = (11) + (H) + (H). (37) 

Next we find the current and the density. They are given by 

-V(^-K-^)^^-^^^^^' ^'^'^ 



27r)'^ V dsks J dfi 

where is the density of states. We find two transport equations for current and charge density or chemical potentials 
by integrating (|3^) and also ( p7\ ) multiplied by Vx with respect to k, and we arrive at Eq. (1), where 

(39a) 
(39b) 

\'tr,s ■'st,tr/ 

^3^,, = !^If!:£. (39c) 
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^V^ti-.s Tsftr 



In (Ua) we have introduced the drag conductivity defined in Eq. (0). In deriving the drag term, we have made use of 
the result obtained for usual Coulomb drag in e.g. Ref. ( |13|). Furthermore, the local electro chemical potential has 
been defined as 

/is = /is + 60, (40) 

where (j) is the electrical potential. 

V. EVALUATION OF THE SPIN DRAG RESISTIVITY 
A. One dimension 

The polarization function is in one dimension at small temperatures, where we can perform an w-expansion, given 

by 

Inixs 9,t^ ~tj— 3— . (41) 

Inserting this into the formula for a^ , performing the ^-integration for the case of a non- magnetic conductor, 
Efcs = £fc, and using that [/o(e)]^ — {&kT)^^5{e — Ef), we find 

The spin-drag resistance is thus proportional to temperature and dependent on the Coulomb backscattering matrix 
element. Clearly, this contribution can be very large at finite temperatures. However, in strictly one dimensions, 
where Fermi liquid theory is not expected to apply, the Boltzmann equation is not a correct starting point and one 
should be somewhat careful about drawing firm conclusions from this. Nevertheless, this Fermi's Golden rule result 
is indicative of e-e interactions being very important for spin transport in one dimension. 

B. Two dimensions 

For the two dimensional case we start by deriving a low temperature result and go on to compare it with a full 
numerical integration of ajj . At small uj and q the imaginary part of the polarization function is given by 

2 
ui 

Imxsiq, ^) « uj 3 , (43) 

2nn q-^kp 



and the screened static Coulomb interaction is 



e2 



Uiq) = t: 7 7- (44) 

2Treoer {q + qxp) ^ ' 



In this approximation, the q-integral becomes 

dQi'^r— — ^' (^^) 

[q + qTFr 

which is clearly not convergent and therefore we set the upper limit to be 2kp , because Im x is zero for momentum 
exchange larger than 2kp. With these inputs, we arrive at the following approximate expression 

where 7 = 2kp /qTF and qrp — TOe^/27reo£r^ is the two dimensional inverse Thomas-Fermi screening length. Typi- 
cally 7 is of order 1. This means that spin drag resistivity can be equal to a fraction of the quantum resistance and 
should therefore indeed be measurable for standard high mobility quantum wells. 

We have also integrated the spin drag formula numerically, see Ref. flj] for details. The result is shown in Fig. n^ 
for realistic numbers for a two dimensional GaAs electron gas. The integration is done using the full dynamically 
screened interaction for a quantum well with finite thickness. The approximate formula (Eq) is seen to overestimate 
the spin drag effect slightly. 
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FIG. 1. The spin-drag resistivity in the two-dimensional case as a function of temperature. The full line is the numerical 
integration of Eq. (Gl) for a two-dimensional quantum well of thickness 10 nm, electron density 2x 10^5 m~^. We have used 
typical parameters for GasAs based heterostructures. The right most thin line is the approximate expression in Eq. (jM), while 
the left thin line is the result of a integrating Eq. (0), but using the T = expression for x('?,'^)- 

VI. CONCLUSIONS 

We have derived a set of transport equations for spin polarized drag which incorporate e-e scattering. This has 
been done within the framework of the Bohzmann equation. First we showed that the in the absence of spin-flip 
scattering and at low temperatures the exact solution of the Boltzmann equation corresponds to two shifted Fermi- 
Dirac distribution functions. Furthermore, if the interaction is weak one can use perturbation theory and arrive at 
the same conclusion following the lines of argument from usual Coulomb drag. Having observed that the shifted 
Fermi-Dirac distribution is correct at low temperatures or weak e-e scattering, we go on to the general case which is 
solved approximately by using the SFD as an ansatz, which allows for a solution of coupled Boltzmann equations. 

The main conclusion from this is that e-e interaction introduces a spin drag term, which tend to drag the spin 
currents to be equal. There a two such mechanisms namely e-e interactions, which is temperature dependent, and 
angular dependent elastic spin flip scattering, which is temperature independent. Therefore if a spin polarized current 
is driven through the system, the spin drag will give rise to an additional resistivity. This resistivity increases with 
temperature. We have solved for the spin drag resistivity numerically in two dimensions which shows that it can 
become considerable and even exceed the ordinary impurity scattering induced resistivity. The spin drag should thus 
be measurable in for example a structure combining a two dimensional electron gas with ferromagnetic materials or for 
one-dimensional systems e.g. made by nanotechnology in semiconductors or by contacting nanotubes to ferromagnetic 
contacts. 
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